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ABSTRACT 

m _ 

We present an analysis of the production and two-photon decay of the lightest 
CP-even Higgs boson of the Minimal Supersymmetric Standard Model (MSSM) at the 
. Large Hadron Collider (LHC). A rather general model is considered, without super- 

gravity constraints. All parameters of the model are taken into account, we especially 
study the dependence of the cross section on the squark masses, on the bilinear pa- 
rameter fj, and the trilinear supersymmetry breaking parameter A. Non-zero values 
Oh, of these parameters lead to significant mixing in the squark sector, and, thus, affect 

the masses of Higgs bosons through radiative corrections, as well as their couplings to 
squarks. The cross section times the two-photon branching ratio of h° is of the order 
of 15-25 fb in much of the parameter space that remains after imposing the present 
, experimental constraints on the parameters. 



1. Introduction 

The most important production mechanism for the neutral SUSY Higgs bosonJH 
at the Large Hadron Collider (LHC) is the gluon fusion mechanism, pp — ► gg — > h°, 
H°, A ,! and the Higgs radiation off top and bottom quarks.! Except for the small 
range in the parameter space where the heavy neutral Higgs H° decays into a pair 
of Z bosons, the rare 77 decay mode, apart from rr decays, is a promising mode 
to detect the neutral Higgs particles, since b quarks are hard to separate from the 
QCD background. It has been pointed out that the lightest Higgs could be detected 
in this mode for sufficiently large values of the mass of the pseudoscalar Higgs boson 
rriA > mz- 
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Here we present results of a recent study-3 of the hadronic production and subse- 
quent two-photon decay of the lightest CP-even Higgs boson (h°) of the MSSM, which 
is valid for the LHC energy of \/s = 14 TeV, focussing on the case of intermediate- 
mass squarks. The case of heavier squarks has been discussed elsewhere. A related 
study has been presented by Kane et al.H They consider a model where parameters are 
related by supergravity, but otherwise chosen randomly within their allowed ranges. 
As mentioned, the gluon fusion mechanism is the dominant production mechanism 
of SUSY Higgs bosons in high-energy pp collisions throughout the entire Higgs mass 
range. We study the cross section for the production of the h° and its decay, taking 
into account all the parameters of the Supersymmetric Standard Model. In particular, 
we discuss the dependence on the squark masses, and take into account the mixing 
in the squark sector, the chiral mixing. This also affects the Higgs boson masses 
through appreciable radiative corrections, and was previously shown to lead to large 
corrections to the rates.i 

In the calculation of the production of the Higgs through gluon-gluon fusion, 
we include in the triangle graph all the squarks, as well as b and t quarks, the lightest 
quarks having a negligible coupling to the h . On the other hand, in the calculation 
of decay of the Higgs to two photons, we include in addition to the above, all the 
sleptons, W ± , charginos and the charged Higgs boson. 

An important role is played in our analysis by the bilinear Higgs coupling /x, 
which occurs in the Lagrangian through the term 



where H\ and H 2 are the Higgs superfields with hypercharge —1 and +1, respectively. 
Furthermore, the Minimal Supersymmetric Model contains several soft supersym- 
metry-breaking terms. We write the relevant soft terms in the Lagrangian as follows^ 



(see also ref. [9]) with subscripts u (or U) and d (or D) referring to up and down-type 
quarks. The Higgs production cross section and the two-photon decay rate depend 
significantly on several of these parameters. 

The Higgs production cross section and the two-photon decay rate depend 
on the gaugino and squark masses, the latter being determined by, apart from the 
soft-supersymmetry breaking trilinear coefficients (A u , An) and the Higgsino mixing 
parameter /x, the soft supersymmetry-breaking masses, denoted in eq. @ by Mu, fhu 
and rfio, respectively. For simplicity, we shall consider the situation where Mu = 



-fiH^eH 2 +h.c., 



(1) 




-MuQ ] Q - mlu R] u R - m 2 D d R] d R - MI x E\E x - M^H^ 

+— IaA + A A j + — {A fe A fc + A k A k } , 
2 2 , , 



(2) 
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fhu = ffi£) = m, with m chosen to be 150 GeV for the first two generations^ and 
varied over the values 150, 250, 500 and 1000 GeV for the third generation. As 
discussed in ref. [6], it suffices to consider A positive and vary the sign of \x. 

In Sec. 2, we discuss the implications of the nonzero values of A and \x on the 
Higgs masses, together with the constraints related to the other relevant masses. We 
then go on to study the cross sections and decay rates for the lighter CP-even Higgs 
boson in Sec. 3. 



2. Constraints on the Parameter Space 

In this section we describe in detail the parameter space relevant for the pro- 
duction and decay of the lightest Higgs boson at LHC, and the theoretical and ex- 
perimental constraints on it before presenting cross sections and decay rates. 

At the tree level, the masses of the CP-even neutral Higgs bosons are given 
by (m h o < m H o )0 



1 

m H°,h° — Q 



m A + m l ^ V ( m A + m l) 2 — 4m|m^ cos 2 2/3 



(3) 



which are controlled by two parameters, rriA (the mass of the CP-odd Higgs bo- 
son) and tan/3 (= t^/fi, where v-i and v\ are the vacuum expectation values of the 
two Higgs doublets). Indeed, the entire Higgs sector at the tree level can be de- 
scribed in terms of these two parameters alone. At the tree level, the mass of the 
lightest Higgs boson (m/jo) is bounded by m^. There are, however, substantial ra- 
diative correctionsOl to the CP-even neutral Higgs masses. The radiative corrections 
are, in general, positive, and they shift the mass of the lightest neutral Higgs boson 
upwards.!! More recent radiative corrections to the Higgs sectoral which are valid 
when the squark masses are of the same order of magnitude, have not been taken into 
account in our studyl As long as the "loop particles" are far from threshold for real 
production, the cross section does not depend very strongly on the exact value of the 
Higgs mass. 

In order to simplify the calculations, we shall assume that all the trilinear 
couplings are equal, 

A u = A d = A. (4) 

Furthermore, we shall take the top-quark mass to be 175 in our numerical 

calculations. The parameters that enter the neutral CP-even Higgs mass matrix are 
varied in the following ranges: 

50 GeV < rriA < 1000 GeV, 1.1 < tan (3 < 50.0, 

50 GeV < < 1000 GeV, < A < 1000 GeV. (5) 



''increasing m to 500 GeV for the first two generations leads to an increase in the cross section times the 
two-photon decay rate of about 5%. 
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These values cover essentially the entire physically interesting range of parameters in 
the MSSM. However, not all of the above parameter values are allowed because of 
the experimental constraints on the squark, chargino and h° masses. For low values 
of rh, the lightest squark tends to be too light (below the most rigorous experimental 
bound of ~ 45-48 GeV0) , or even unphysical (mass squared negative). The excluded 





Figure 1. Regions in the /i-tan/3 plane which are ruled out by too light chargino (x^) and squark 
masses. The gaugino mass scale is M 2 = 200 GeV and tua = 200 GeV. The solid (dashed) contours 
for small refer to the chargino mass rn x ± = 68 (90) GeV. Two values of rh are considered, 
left: rh = 150 GeV, right: rh = 250 GeV. For each value ofrh, two values of the trilinear mixing 
parameter are considered. For rh — 150 GeV, the squark masses are too light or unphysical in 
much of the /i-tan/3 plane. The hyperbola-like contours give regions that are excluded by the 
lightest b squark being below 45 GeV (solid) or 90 GeV (dashed). The more straight contours at 
large fi and small tan/3 similarly indicate regions that are excluded by the lightest t squark. 

region of the parameter space is shown in the left part of fig. 1 for rh = 150 GeV, 
M2 = 200 GeV, rriA = 200 GeV and for two values of the trilinear coupling A. The 
allowed region in the /i — tan/3 plane decreases with increasing A, but the dependence 
on M 2 and in this region is rather weak. In order to have acceptable 6-squark 
masses, /x and tan/3 must lie inside of the hyperbola-shaped curves. Similarly, in 
order to have acceptable t-squarks, the corners at large |/x| and small tan/3 must be 
excluded. 



4 



The chargino masses are, at the tree level, given by the expression 
1 



m 2 ± = -{M 2 + fi 2 ) + m 2 w 



± 



n -i 1/2 

-{Ml - fi 2 ) 2 + r/4 cos 2 2(3 + m 2 w {Ml + ii 2 + 2/iM 2 sin 2(3) . (6) 

When n = 0, we see that, for tan/? 3> 1, the lightest chargino becomes massless. 
Actually, small values of fi are unacceptable for all values of tan/?. The lowest ac- 
ceptable value for \fi\ will depend on tan/?, but that dependence is rather weak. The 
excluded region due to the chargino being too light, increases with decreasing values 
of M 2 . We note that the radiative corrections to the chargino masses are small for 
most of the parameter spaced In fig. 1 we show the contours in the /z-tan (3 plane 
outside of which the chargino has an acceptable mass (> 68 GeV)lll. By the time 
the LHC starts operating, one would have searched for charginos with masses up to 
90 GeV at LEP2. Contours relevant for LEP2 are also shown. 

For m = 250 GeV (right part of fig. 1), the region excluded due to low or 
unphysical squark masses is much reduced, and at m = 500 GeV it is practically 
absent.! However, for larger values of m the experimental constraints on the h° masslll 
rule out some corners of the /i — tan (3 plane.H The extent of these forbidden regions 
in the parameter space grow rapidly as itla decreases below (9(150 GeV). They also 
increase with increasing values of A. 

As discussed above, the mass of the lighter CP-even Higgs boson h° will depend 
significantly on fi, tan/3, A and in, through the radiative corrections. For m = 
250 GeV, and two values each of m A (100 and 200 GeV) and A (0 and 1000 GeV), 
the dependence on /i and tan/3 is displayed in fig. 2. At large and large tan/3, the 
radiative corrections are large and negative, driving the value of m^o well below the 
tree-level value. (Those regions practically coincide with those where the b squarks 
are very light or unphysical.) 

The charged Higgs boson mass is given by 

m 2 H ± = m 2 w + m 2 A + A, (7) 

where A arises due to radiative corrections and is a complicated function of the 
parameters of the model. 

The radiative corrections to the charged Higgs mass are, in general, not as 
large as in the case of neutral Higgs bosons. In certain regions of parameter space 
the radiative corrections can, however, be large. This is the case when the trilinear 
mixing parameter A is large, itla is small, and when furthermore tan (3 is large. We 
shall include the effects of non-zero A and ix in the calculation of the charged Higgs 
mass. The present experimental lower bound of 40-45 GeV^ on the charged Higgs is 
not restrictive, but presumably by the time the LHC starts operating, one will have 
searched for charged Higgs bosons at LEP2 with mass up to around 90 GeV. Even 
this bound does not appreciably restrict the parameter space. 
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Figure 2. Mass of the lightest CP-even Higgs boson vs. fi and tan/3, for M2 — 200 GeV and 
rh — 250 GeV. Two values of m A and two values of A are considered: a) m A = 100 GeV, 
A = 0, b) m A = 200 GeV, A = GeV, c) m A = 100 GeV, A = 1000 GeV, d) m A = 200 GeV, 
A = 1000 GeV. 

The neutralino mass matrix depends on the four parameters M 2 , Mi, fi and 
tan/?. However, one may reduce the number of parameters by assuming that the 
MSSM is embedded in a grand unified theory so that the SUSY-breaking gaugino 
masses are equal t o a common mass at the grand unified scale. At the electroweak 
scale, we then have&H M\ = (5/3) tan 2 #w M%. We shall assume this relation through- 
out in what follows. The neutralino masses enter the calculation through the total 
width of the Higgs boson. We here present numerical results for the gaugino mass 
parameter being M2 = 200 GeV. A wider range of values is considered in ref. [6]. 
The experimental constraints on the lightest neutralino mass rule out regions of the 
parameter space similar to those ruled out by the charginos.il 
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3. Cross Section and Two-Photon Decay 

The cross section for pp — * h°X, is calculated from the triangle diagram con- 
voluted with the gluon distribution functions, 



„ = V2,G F Of I E/,Wf /'° g,V " /m " 0> d!/G (^e») G(^ 



with contributions from various diagrams (k). For the standard case of a top-quark 
loop, 



T 



^l-(r-l) 



arcsm 



1 



(9) 



and r = (2m t /mho) 2 > 1. 

For M 2 = 200 GeV, m = 500 GeV, and fi = 500 GeV, plots are given in 
refs. [6,7]. The following features are noteworthy: 

• The cross section decreases appreciably for large values of A. This is mainly 
due to an increase in the h° mass. 

• The cross section increases sharply for small values of tan (3, and also at small 
rriA- The increase at small tan/3 is caused by the h° becoming light. At small 
values of tua and large A, the couplings of h to b quarks and r leptons become 
large, making the cross section very large in this region. 




Figure 3. Cross section for pp — ► h°X — > jjX as a function of wia and tan (3 for = 200 GeV, 
two values ofm: left: m = 500 GeV, right: m = 1000 GeV. For both cases, /i = 500 GeV, and 
A = Q. 

The two-photon decay rate is found! to increase sharply at large values of 
A, but this does not result in a correspondingly larger rate for the process 

pp -> h°X -> 77X, (10) 
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since the production cross section also decreases. In fig. 3 we show the cross section 
for the process (|10"D for the case of heavy squarks. A characteristic feature of the cross 
section is that it is small at moderate values of ttia, and then increases steadily with 
increasing tua, reaching asymptotically a plateau. This behaviour is caused by the 
contribution of the W to the triangle graph for h° — *> 77. The h°WW coupling is 
proportional to sin(/3 — a), where 

cos^-a) = m rj m }-'% (11) 
m A o{m z H o - m^o) 

For large m^, at fixed (3, all Higgs masses, except m h o, become large, so that h° 
decouples. For large m^, we actually have sin(/5 — a) — > 1, which is why the cross 
section increases and reaches a plateau for large m^. 




Figure 4. Cross section for pp — > h°X — > jjX as a function of uia and tan (3 for M2 — 200 GeV, 
rh = 250 GeV, n = 500 GeV, and two values of A: left: A = 0, right: A = 200 GeV. 

The situation changes dramatically when we consider lighter squarks. In fig. 4 
we show the corresponding cross section for the case of m — 250 GeV. The main 
difference is that for large values of tan/? the product of the cross section and the 
two-photon decay rate practically vanishes. This is due to the fact that b squarks 
become light in this limit, the h° can decay to b squarks, and the two-photon decay 
rate becomes too small. 

For A = 200 GeV, there is a band of larger values at tan (3 ~ 26. This region 
is rather "turbulent" (for these values of rh and A): As mentioned above, the two- 
photon decay rate vanishes due to decays to light b squarks, but at the same time the 
light Higgs tends to make the cross section big. The competition between these two 
effects is responsible for the turbulent behaviour seen here. 
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For some parameters, the latter effect may dominate, and we get spikes or 
bands in the product of the cross section and the two-photon decay rate. (Since 
a small Higgs mass is due to large radiative corrections — which are not accurately 
known — it is not clear how physical these bands or spikes are.) 

The /i-dependence of the cross section is for the case of M 2 = 200 GeV and 
m = 500 GeV discussed in refs. [6,7]. Here, for m = 250 GeV, we show in fig. 5 the 
cross sections corresponding to those of fig. 4, but for \i = —500 GeV. 




Figure 5. Cross section for pp — > h X — > 77 V as a function of uia and tan (3 for M2 = 200 GeV, 
m = 250 GeV, ,u = -500 GeV, and two values of A: left: A = Q, right: A = 200 GeV. 



4. Summary and concluding remarks 

We have discussed the cross section for the production of the lightest CP-even 
Higgs boson at the LHC, in conjunction with its decay to two photons. Where the 
parameters lead to a physically acceptable phenomenology, the cross section multi- 
plied by the two-photon branching ratio for the lighter CP-even Higgs boson is of the 
order of 15-25 fb. 

Similar results have been presented in ref . [8] . Within the context of a SUGRA 
GUT model, these authors consider basically a random sample of parameters com- 
patible with experimental and theoretical constraints. The cross sections obtained in 
ref. [8] appear to be somewhat higher than those of ref. [6]. 

It should be noted that in regions where the Higgs cross section times the 
two-photon decay rate is small, typically the lightest b squark is light. Thus, as 
"compensation" , one should be able to observe b squarks. 

The recent Fermilab data on large E t inclusive jet cross sectionsil suggest that 
the gluon distribution functions are larger at high x.c3 We have checked whether these 
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distributions (specifically, CTEQ4H J) lead to a higher rate for Higgs production. This 
turns out not to be the case. The enhancement of the gluon distribution function is 
in a range of x where its magnitude is simply too small, anyway. 

These calculations do not take into account QCD corrections. Such corrections 
have been evaluated for the quark-loop contribution, and lead to enhancements of the 
cross section of about 50%.E-3 For the squark loops QCD corrections have recently been 
studied using low-energy theorems.!^ It is concluded that the additional contributions 
lead to the same QCD corrections as for the top and bottom quark loops. 

It is a pleasure to thank the Organizers of the Minsk Workshop, in particular 
Professor L. Tomil'chik, for creating a very stimulating and pleasant atmosphere 
during the meeting. This research has been supported by the Research Council of 
Norway. The work of PNP is supported by the Department of Science and Technology 
under project No. SP/S2/K-17/94. 
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